A model for single electron decays from a strongly isolated quantum dot 
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Recent measurements of electron escape from a non-equilibrium charged quantum dot are inter- 
preted within a 2D separable model. The confining potential is derived from 3D self-consistent 
Poisson- Thomas-Fermi calculations. It is found that the sequence of decay lifetimes provides a 
sensitive test of the confining potential and its dependence on electron occupation. 

73.23.Hk, 73.20.Dx, 73.40Gk, 85.30.Vw, 71.10.-w 

I. INTRODUCTION 
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In a recent experiment |jj , a strongly isolated quantum 
dot was charged with excess electrons, and their sequen- 
tial escapes were recorded over a one hour time period. 
This was repeated 150 times to obtain a statistical distri- 
bution of decay times. The dot is formed in an electron 
gas located at a depth of 70 nm in a GaAs — AlGaAs het- 
erostructure. Its shape is defined by electrostatic confine- 
ment using a set of gates, as sketched in the insert to Fig. 
[I] . The gate voltages were ramped up quickly, so that the 
dot retained a sizeable number of excess electrons when 
it was well isolated from the surrounding electron gas. 
The observations correspond to sequential tunnelling of 
(seven) electrons from the dot to the surroundings. The 
lifetimes extracted from the escape times distribution |l]] 
are shown in Fig. |l|. A striking quasi-linear dependence 
of the logarithm of the lifetime on electron number is 
apparent. 

Sequential decays have been known and studied for 
over a century in the context of nuclear physics. The 
combined instances of alpha and beta decays from the 
heaviest elements are responsible for most natural ra- 
dioactivity. The description of alpha decay in terms of 
tunnelling of alpha particles through a confining potential 
dates back to the 1920's (Gamow ||, Condon and Gurney 
||). Although the basic nature of the decay as a barrier 
penetration is well understood, accurate predictions for 
radioactive lifetimes are difficult because the process by 
which the escaping alpha-particle is preformed within the 
nucleus requires an understanding of four-body correla- 
tions. As a result, it is impossible to deduce accurate 
information on the barrier shape. Nevertheless, alpha 
decays have provided useful information on nuclear radii 
and the range and gross features of the nuclear interac- 
tionfll. 
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FIG. 1. Experimental lifetimes (in seconds) extracted 
from the decay sequences, as reported ref. [1]. The insert 
shows the gate arrangement which defines the dot. 

It has become commonplace to say that a quantum 
dot is an artificial atom, but in fact the self-consistent 
potential confining electrons in a large dot has more in 
common with the mean field potential in a heavy nu- 
cleus: flat in the interior, with abrupt walls. An artificial 
nucleus is a more apt description, as will become clear 
in this paper. Indeed, the detection of sequential decays 
from an isolated quantum dot is a more favourable sit- 
uation for study of the decay process, as the question 
of preforming the electron does not arise. Hence, we can 
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more confidently test our knowledge of the confining bar- 
riers for electrons, as well as the profile, and dependence 
on occupation number, of the dot potential. We will an- 
alyze these aspects in this work, and show that these 
measurements of the lifetimes of "radioactive quantum 
dots" introduce new constraints on our ability to model 
their structure. 

The present experiment has another significant advan- 
tage over nuclear decays: instead of counting incoherent 
decays from a large sample of identical nuclei, here a 
single dot is involved, and the correlation between con- 
secutive events can be analyzed. In addition, it should 
be possible to design the shape, density and excitation 
energy of the dot within rather broad margins, so that 
future experiments on mesoscopic systems will be much 
more flexible than those in nuclear systems, where only 
those nuclei existing in nature, or created in sufficient 
numbers, can be studied. Thus, the study of electron 
decays from a quantum dot has the potential to reveal 
new features of the tunnelling process. This is a topic of 
currently renewed interest: see for example van Dijk and 
Nogami || . The type of simple model developed in this 
paper can be of great utility in such future studies. 

In this work we will describe the decay process using 
analytic models which incorporate characteristics of the 
confinement potential extracted from realistic numerical 
simulations. As the dot contains about 300 electrons, 
Poisson-Thomas-Fermi calculations should be adequate 
to describe the electron density and the confining poten- 
tial of the dot. With these in hand we have developed 
accurate analytic approximations for the confining poten- 
tial that allow us to construct an envelope approximation 
wavefunction for the electrons in the dot, and to compute 
the electron lifetimes from a fully quantal expression for 
the transmission amplitude across the barrier. 

Previous works which model a quantum dot have been 
concerned with the wave functions of confined states in 
the dot, the electron density distribution and the shape 
of the confining potential. For such purposes, only the 
inside of the barrier matters. It is when one looks at the 
escape of electrons from the dot that the barrier height, 
its width and shape become important; these are the 
new features explored in this paper. In section II we 
describe the development of our model, while in Section 
III we discuss the results for the sequence of lifetimes 
and compare them with experiment. Some details are 
relegated to two appendices. 

II. MODELLING OF ISOLATED DOT DECAYS 

A. Framework 

The Poisson-Thomas-Fermi modelling is described in 
more detail elsewhere ||, so here we list only the main 
steps: 



i) first, Poisson-Schrodinger (PS) and Poisson-Thomas- 
Fermi (PTF) simulations as described in jjj are per- 
formed for the ungated heterostructure. Our inputs for 
the PS simulation are the thickness and composition of 
each layer in the heterostructure, and the dopant con- 
centration in the donor layer. From these we predict 
the density of the 2DEG. The only adjustable param- 
eter is the donor ionization energy which is set to be 
e$j = 0.12 eV, in order to reproduce the measured 2DEG 
density, n e — 2.74 10 11 cmT 2 . For the simpler Poisson- 
Thomas-Fermi scheme we employ a common relative per- 
mitivity e r = 12.2 for all layers of the heterostructure, 
which, combined with the parameters already used for 
the PS simulation, also reproduces the experimental n e . 
After this "fitting" the model has no other free parame- 
ters. 

ii) For the gated structure we use the gate layout and 
voltages of the experiment. To solve the Poisson equa- 
tion for the gated heterostructure one has to impose as a 
boundary condition the value of the electrostatic poten- 
tial on the exposed surface of the heterostructure, and 
on the gates. We assume Fermi level pinning and choose 
the energy of the surface states as the zero of the energy 
scale. In this convention, the conduction band edge is 
set at eV s — 0.67 eV on the exposed surface. Under each 
gate the conduction band is set at eV ms + eV g , where V g 
is the gate voltage and the metal semiconductor contact 
potential, eV ms , is taken as 0.81 eV The electrostatic 
potential due to the gates is then computed using semi- 
analytic expressions based on the work of Davies et al. 
]l0| and pT] . Added to this are: a) the Coulomb po- 
tential (direct term) between the electrons, and a mirror 
term which imposes the boundary conditions at the sur- 
face, and b) the contribution from the fully ionized donor 
layer and its mirror term (see Sect. IIA of |sj for details 
of a similar example.) We neglect exchange and correla- 
tion effects, which are small. 

Hi) The connection between the confining potential de- 
fined by the conduction band edge and the electron den- 
sity is completed by using the Thomas-Fermi approxima- 
tion at zero temperature: 

Pe(r) = ^ (^3-(^ " ^«)) 3/ " (1) 

The PTF iteration is performed starting from the un- 
gated heterostructure densities as trial values. 

B. Equilibrium dot 

As a first step, we examine the dot in its final state 
after all the excess electrons have escaped. This corre- 
sponds to a PTF simulation with the same Fermi level, 
EF,dot = 0, for the electrons in the dot and in the 2DEG 
outside the barriers. The gate voltages are taken from 
ref. f§ as V PL = -0.40 V, V C1 = V C 2 = -0.44 V and 
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V H = -0.7 V. The predicted PTF 3D electron distribu- 
tion p e (x,y,z) is more conveniently visualized in terms 
of a projected 2D density: 

poo 

n e (x,y) = p e (x,y,z)dz, (2) 

J Zj 

where Zj is the junction plane. The n e (x,y) distribu- 
tion, shown in Fig. ^, has an approximately rectangular 
boundary, and its maximum value is close to the 2DEG 
density of the ungated heterostructure. In this calcula- 
tion the dot contains 286 electrons. 




Although the computed V(x, y, z) is not separable, pre- 
vious experience with Poisson-Schrodinger simulations of 
wires (l2) and circular dots, has shown us that the 
factorization ansatz leads to very good approximations 
when the z degree of freedom is integrated out as in eq. 
||. This prescription to construct the 2D potential avoids 
the type of ad hoc assumptions often made. 

In Fig. H we show the UpTF(x,y) corresponding to 
the equilibrium dot of Fig. |[ As expected from the 
gate layout shown in the inset to Fig. ffl, it has two very 
high barriers running parallel to the x axis, one centered 
at y = and the other that begins with a steep rise at 
y ~ 400 nm (and shows clearly the mark of the three- 
fingered gate layout labelled CI, C2, and PLunger in 
Fig. |l|). Tunnelling across these barriers is negligible. 
In addition there is a symmetric pair of barriers running 
parallel to the y axis, with maxima at x ~ ±238 nm 
through which the electrons do tunnel. In the interior, 
the potential is practically constant. Although these x- 
barriers have somewhat increasing height with increasing 
y, the rectangular shape of the potential suggests using 
a separable approximation in cartesian coordinates: 

U PTF (x,y) ~U s {x,y) = U(x)+W(y) (5) 



FIG. 2. The two dimensional PTF density, n e (x,y), for a 
dot in equilibrium with the surrounding 2 DEG. 



C. Dot with excess electrons 

To study these configurations we set the Fermi level in- 
side the dot, Ep^oti higher than its value outside the bar- 
riers, Ep,2DEG — 0- We can do so because the dot is well 
pinched off from the surrounding electron gas. We ran 
PTF simulations with equally spaced values for Ep^ot 
running from to 17.5 meV in steps of 2.5 meV. The 
occupation Q of the dot increases linearly with Ep^ot 
at the rate 2.75 electrons per meV, giving occupations 
286 < Q < 334. 

The simulations also produce the confining potential 
for the electrons in the dot, eV(x, y, z). To reduce this to 
a two-dimensional function, U(x,y), we take a weighted 
average over the density profile in the z direction: 



where 



UpTF{x,y) 



P(z) 



f™eV(x,y,z) P(z) dz 



p e (x,y,z) dx dy 



(3) 



(4) 



Here the domain of integration D, is a rectangle in the xy 
plane which extends a short distance into the surround- 
ing electron gas, (from (xi,yi) = (—510 nm, —255 nm) to 
(x r ,y r ) = (510 nm, 255 nm).) This includes an area out- 
side the dot where the 2DEG is still depleted by the gates. 
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FIG. 3. Two dimensional confining potential, UpTF(x,y) 
for the dot of Fig. 1. 
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FIG. 4. Continuous lines: sections at y = 200 nm of the 
Uptf{x,v) corresponding to Ep^ot = 0.0,0.005,0.010 and 
0.015 eV. Dashed lines: analytic parametrization for U(x) 
as described in the text (the latter shown only for x > for 
clarity.) 

We will interpret the experimental decay data using 
this separability ansatz. For the W(y) barriers, which 
are basically impenetrable, we use two simple models de- 
scribed below. As a guide to a realistic choice for the 
x-dependent term we examine in Fig. ^| the profiles of 
Uptf{x, y) at a fixed value of y = 200 nm in the middle 
of the dot. The profiles shown cover a range of occu- 
pations of up to forty excess electrons. In this range, 
the potential at the dot center increases linearly with Q, 
according to 



Uq = 0.347Q - 118.4 meV 



(0) 



At large distances outside the dot, U^ = —18.8 meV is 
constant. Similarly, the location of the barrier maximum 
and its height can be parametrized as: 



Q-286 

xi, = Z6o nm 

16 

U b = 0.U7Q - 13.4 meV 



(7) 



Note that dUb/dQ » l/3dUo/dQ reflects the decrease 
of the screened Coulomb repulsion away from the cen- 
ter of the dot. Furthermore, we have found that the 
x-dependence can be very well reproduced (see Fig. 
using the following analytic model: 

U(x) = U b + U M F(x), x>0, 

= U{—x) , x < 0, where 



sinh 2 



Umf = U c 



I x-x b \ 
\ w b J 



cosh 



(8) 



This potential form has the great advantage that the 
transmission coefficient for Umf 1S known analytically 
p3[ . Umf is an asymmetric barrier which takes one value 
for x « xi, and another value for x >> Xb. 



U M F(xb) = 

Umf(oo) = lim Umf(x) = U c e 2 ^ 

x — >oo 

Umf(-oo) = lim Umf(x) = U c e~ 2fl 



(9) 



The parameters Ub, U c , fj,, Xb, Wb allow one to fit the bar- 
rier height, the potential floors inside and outside the 
dot, the barrier spacing and the barrier width. Since the 
barriers are spread quite far apart, in practice Xb >> Wb, 
so Umf(x = 0) w Umf(~ °°)- I n this case, 



U = 17(0) « U b + U c e- 2 » 
Uoo ee lim U{x) = U b + U c e 2f *. 



(10) 



Then we can solve for 



1, (Ub-Uoo 



and 



U c = -{U b ~ U Q )e 



2/J 



(11) 



To determine the parameters appearing in eq. [8|, we take 
the values of the PTF potential at the origin, Uq, well 
beyond the barrier, Uoo, and the value U Xb at the barrier 
maximum x = Xb, and then plot U(x) to find the best Wb, 
which turned out to be 48 nm. This gives a convenient 
analytic form for the confining potential, motivated by 
PTF, whose transmission coefficient is: 



2 sinh(7r£; + ) sinh(7rfc_ ) 
cosh(7r(fc + + &_)) + cosh(7r/3) 



(12) 



where 



k -/+ = 



2m* 



(E - U 0/oo )w 2 



and 



2m* 



(2U b - 2U C -U Q - Uoo)v% - 1 



(13) 



Barrier shape W(y): In Fig. [5] we examine a section 
of Uptf(x = Q,y) through the center of the dot. We 
use two approximate models, the simplest one being an 
infinite square well, of width w y ~ 350 nm. The slightly 
fancier one is a truncated harmonic oscillator: 



W tho {v) 







(flat bottom) 



-0.13 + -k y (y - y ) 2 (walls) 



(14) 



with yo = 238 nm and k y = 7.35 ■ 10~ 6 nm~ 2 . As can 
be seen in Fig. || this parametrization (plus the constant 
term Uq) reproduces the main features of the x = sec- 
tions of the PTF potentials. 

By combining eqs. 6] to |l4 we determine a separable 
analytic potential model for the dot containing a desired 
number Q of electrons. This removes the necessity of re- 
peatedly solving the PTF equations for the self-consistent 
field, while studying the decay process. 
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FIG. 5. Continuous lines: sections at x = nm of the 
Uptf(x,v) for the same EF,dot as in Fig. 3. Dashed lines: 
analytic parametrization of W t ho(y) + Uo as described in the 
text. 



D. Quasibound states of the dot 

We construct the electron wave functions inside the 
dot in the envelope function approximation, using our 
parametrized potential, U s (x,y). The single electron en- 
ergies are 



En x ,n v E nx + E n 



and the electron wavefunctions factorize 

^n x ,n y {x,y) = (l>n x (xi)lp nv (y) ■ (16) 

The factors satisfy ID Schrodinger equations: 

:4>L{ X ) + U { X )^A X ) = E n x </>n x (x) 



2m 



2m 



<M + W(y)i> n «(y) = E n J ny (y) (17) 



The second equation is for a confined wavefunction, eas- 
ily solved by standard numerical methods. We label the 
solutions by the number of loops, n y , of the eigenfunc- 
tion. For example, taking W(y) to have hard walls, the 
energy is 



*2 / \ 2 

n I riyir 



2m* V w 



J y 



(18) 



For the truncated harmonic oscillator shape there is no 
similar analytic expression, but the dependence on n y is 
similar. 

The x-dependent equation describes ID electrons con- 
fined in the dot by the "leaky barriers" . Weakly qua- 
sibound state solutions were computed using methods 
described in However, for levels corresponding to 

the long tunnelling lifetimes observed in the experiment, 
the energies and eigenfunctions can be computed well 



enough by the simpler prescription of setting the electron 
wavefunction to zero at the points ±Xf, inside the barri- 
ers. Furthermore, if only the eigenvalues and lifetimes 
are needed, we have checked that the WKB quantization 
condition is adequate: 



2m* 

—2-{E(n x )-U{x)) dx 



\ 1 ' 



(19) 



In the Appendix we describe the determination of the 
lifetimes, T nx . From here on the energies presented are 
obtained in the WKB approximation. The differences 
from the more accurate predictions using the true quasi- 
bound state energies can scarcely be seen on the scale of 
the graphs. For barrier penetrability we use eq. p"2[ 

We "construct" the desired dot configuration with ex- 
cess electrons by generating a U s (x,y) for the chosen 
value of Q, and filling the levels as follows: a) First we 
list the (E nx , T n ^ ) , in order of increasing n x (and there- 
fore of increasing energy and decreasing lifetime.) This 
list is truncated at an n x — n XtTnax whose lifetime is less 
than 0.01 sec. b) Next we form a list of 2D levels (n x ,n y ) 
by choosing those for which 



E n 



E n „ < E n 



E„ 



(20) 



The levels in this list are occupied in order of increasing 



ene rgy and according to Fermi statistics, see eqs. A5 



(15) A6. We choose the dot Fermi level so that the number of 
electrons is the desired Q. It is supposed that, for the long 
lifetimes observed in the experiment, the electrons have 
time to lose energy by phonon collisions and occupy the 
quasibound states of lowest energy. Then, as described in 
Appendix A, we determine the lifetime for one electron 
to escape from the dot. This involves a weighted average 
of the level lifetimes, according to the occupancy of each 
level at the experimental temperature, T" = 100 mK. 

To produce a sequence of decays for comparison to 
experiment we proceed as follows: i) we start with a 
dot containing a number of electrons, Qo, chosen large 
enough so that the lifetime for one electron to escape 
is smaller than those observed in experiment, ii) We 
redetermine the barrier and dot configuration for Q — 
Qo — 1 electrons, as described in the above paragraph 
and determine again the corresponding lifetime for escape 
of one electron. This process is repeated to generate a 
sequence of decays that covers and extends beyond the 
range of lifetimes measured in experiment. From that 
list we choose as the first observed electron decay that 
corresponding to the Q whose lifetime is the first to be 
larger than to = 25 seconds. 
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FIG. 6. Calculated lifetimes (in seconds) when W(y) is 
either the truncated harmonic oscillator: stars, or a square 
well with w y — 380 nm: + signs . The dotted line is the 
prediction of the two level model, eq. E2. 



III. RESULTS AND DISCUSSION 

In Fig. ||we show results from our model, using param- 
eters chosen as described above, for a range of lifetimes 
extending over three orders of magnitude. The stars cor- 
respond to the truncated harmonic oscillator choice for 
W(y), whereas the +'s are for the square well choice 
(with a value w y — 380 nm chosen to optimize the agree- 
ment with the other prescription in the range of experi- 
mental lifetimes, from 10 to 1000 seconds.) One sees that 
the trends are very similar. For Q in the neighbourhood 
of 304, the predicted decay lifetimes fall in the experi- 
mental range. 

As already mentioned in Section II, our PTF simula- 
tions predict Q — 286 for the dot in equilibrium with the 
surrounding electron gas. This is also what we find with 
this separable model, as the curve of lifetimes shown in 
figure H extrapolates smoothly up to Q = 287, for which 
we predict a lifetime of log 10 r = 5.2, or 44 hours. After 
that, the Fermi level of the electrons inside the dot falls 
below that of the surrounding 2DEG and further decays 
are blocked. It should take almost two days for the dot 
to reach equilibrium with its surroundings. 

Before attempting a more detailed comparison with 
the experimental data it is useful to examine the main 
features in our predicted sequences. First we focus on 
the linear behaviour for values Q < 300. (We have 
found similar behaviour in other ranges of Q when we 
use slightly different sets of parameters.) Such linear de- 
pendence occurs when our model produces a sequence 
of decays dominated by those from a single ID elec- 
tron level; i.e. corresponding to a fixed value of n x . To 
understand why, suppose that at zero temperature and 
for Q electrons, the occupied level with shortest lifetime 
is (n x ,s,n v ), and that {n' x ,n' y } are occupied levels with 
higher energy and longer lifetime (this requires that at 



least n' x < n x ,s for longer lifetime and n y > n y for higher 
total energy.) When one forms the Q — 1 electron config- 
uration according to the rules explained above, one of the 
{n' xl n' y } levels will be empty, whereas the level (n liS , n y ) 
will again be filled. In more physical terms: all the elec- 
trons with energy above that of the level with shortest 
lifetime will lose energy by phonon collisions and fall into 
the leaky level, from which they finally escape. Since the 
lifetime does not depend on n y , all the electrons with 
energy above that of the state {n x ^ s ,n y = 1) will escape 
through the same leaky ID level, n IjS , which remains the 
favoured decay channel as long as it is occupied. There- 
fore the total probability for one electron to escape from 
the occupied states with quantum number n XiS is the 
probability for a single ID electron with energy E nx s , 
multiplied by the number of electrons in occupied states 
with the same quantum number n x s '- q Hx : 



iQ) 



XQ) 



(21) 



and when the occupation g„ x s of the leaky level is con- 
stant, the linear variation of log 10 (T) reflects that of the 
lifetime of the leaky level. This is where the 2D nature 
of the quantum dot asserts its presence, even though the 
decay appears to proceed only in one dimension. 

In Fig. [f] we show the occupations of the two levels 
with the shortest lifetimes. One sees that when Q < 300 
the occupation of the n x = 13 level stays practically con- 
stant and n x ~ 14 level remains empty. For higher values 
of Q both levels contribute significantly to the escape life- 
time. In this situation: 



iQ) = 



1 



q a (Q) gb(Q) 
T a {Q) n(Q) 



(22) 



This is shown as the dotted curve in Fig. ^, and it ac- 
counts very well for the trend of the lifetimes predicted 
by the separable model. 




G 



FIG. 7. Total occupation of levels with n x — 13 (black 
squares) and n x = 14 (open squares). The lines are drawn to 
guide the eye. The truncated harmonic oscillator model was 
used for W(y). 

Our separable model favours the appearance of the 
linear decay sequences, because of the degeneracy in 
lifetime of states with the same n XiS . A non-separable 
model would lift that degeneracy and then the lifetime 
sequences should show a behaviour intermediate between 
the two situations discussed above. In particular, the 
sudden change of slope at Q — 302 in Fig. || would pre- 
sumably spread over a wider range of values of Q. Not 
surprisingly, the predicted lifetimes for the observed de- 
cays depend sensitively on details of the barrier shape. 
Those shown in Fig. |^ correspond to the square well 
choice for W(y), and our standard set of parameters. In 
addition we show how the lifetimes vary when the barrier 
width is changed by amounts ranging from +4% to —3% 
(from left to right). As can be seen, the exact value of 
each decay lifetime depends quite strongly on the bar- 
rier width, as expected for a tunnelling process. But the 
number of slow decays is much more stable: four or five 
in most of the cases shown, and in several cases their life- 
times are quite compatible with the experimental points. 
In particular it is remarkable that a 2% increase in the 
standard barrier width produces a sequence (third line 
from left) in excellent agreement with experiment (dis- 
connected points shifted to extreme left). 

There is a clear distinction between the lifetime trends 
of the thicker and thinner barrier widths. In the latter 
one sees very clearly the transition between escape from 
the n x = 13 and the n x = 14 levels at Q — 302. For 
the thicker barrier widths, escape is dominated by the 
ii x = 14 levels that become progressively more occupied 
above Q — 302. 




Q 

FIG. 8. Slow decays: Calculated lifetime sequences cor- 
responding to variations of the standard barrier width from 
-3% (right) to +4% (left) in steps of 1% . The + signs 
joined by a continuous line (to guide the eye) correspond to 
the prediction for the standard set of parameters. Experimen- 
tal points (left triangles) taken from Fig. 1 with the origin of 
Q shifted arbitrarily. 
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FIG. 9. Same as Fig. 8, but for the fast decays. 

We have explored the dependence of the model pre- 
dictions on changes by similar percentages of the barrier 
heights, potential floor [/ , and the width, w y , of W(y). 
The results are qualitatively similar to those shown above 
for the changes in the barrier width, with a number of 
slow decays ranging from 4 to 6, and in some cases they 
are very similar to the data in Fig. [l]. We therefore 
conclude that our model predictions are quite consistent 
with the experimental trends, although a quantitative 
comparison with the measured lifetimes is hampered by 
the strong sensitivity of tunnelling to any small change 
in the barrier shape. 

Finally we show in Fig. ^ predictions for the fast de- 
cays: their number and location in a graph such as that 
of figure H depends very sensitively on the time (t ) be- 
yond which the experiment measures lifetimes. That is, 
as the dot is isolated there must be a burst of very short 
lived escapes, but after some seconds one reaches the 
stage where separate events can be recorded and the life- 
times deduced. The two dashed horizontal lines in Fig. 
^ correspond to values of log 10 (25 s) and log 10 (35 s). As 
can be seen, when we explore the same range of barrier 
widths as in Fig. |L the number of fast decays above that 
to varies from 3 to 4. The overall trend seems to be con- 
sistent with experiment, in particular if the value of to is 
increased towards the more pessimistic estimate of 35 s. 

IV. SUMMARY AND CONCLUSIONS 

Electron escape from a strongly isolated dot with ex- 
cess electrons has been studied in the framework of the 
self-consistent Poisson-Schrodinger and Poisson-Thomas- 
Fermi approximations. Based on these calculations 
a rectangular separable potential model has been de- 
vised which incorporates the main features of the self- 
consistent field. Rearrangement effects are taken into 
account by recalculating the confining potential U s (x,y) 
after each electron escape. 
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The use of a separable potential introduces certain cor- 
relations in the energy spectrum of the single electron 
orbitals. A more realistic confining potential would have 
a more rounded shape, which would remove the separa- 
bility, and modify those correlations. In the same vein, 
the tunnelling in our simplified model is ID, whereas the 
actual process is 2D. 

We find it quite remarkable that despite all these sim- 
plifications the predictions turn out to be so satisfactory. 
The model therefore may be reliable for extrapolating to 
longer times. For instance we find that the isolated dot 
would hold one excess electron for as long as 44 hours. 
On such a time scale, one could use well isolated dots 
containing a few long lived electrons, to study their en- 
tangled states. This would open an interesting new ap- 
proach to the implementation of quantum computation 
in semiconductor devices. 
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APPENDIX A: LIFETIMES 

We summarize here the expressions relating the life- 
times to the probability of transmission across the bar- 
rier. We follow the standard treatment and definitions 
for alpha decay in nuclear physics, as can be found for 
example in ref. Q. 

Our potential U s (x,y) is separable, and the electron 
can escape only across the barriers in the x-direction. 
Therefore, we have adapted the expressions derived in 
to the ID situation. 

The lifetime r = 1/A is the inverse of the "decay con- 
stant" , defined as the number of "decays" per second per 
parent "dot" . For one dot the electron wavefunction is 
normalized to unity over the volume inside the barriers, 
and A for a given level is just the outgoing flux at large 
distance. 

When the decay probability is small, one can treat the 
electron as confined in the dot. Classically, its trajectory 
will oscillate between the right, x r , and left, xi turning 
points, with a period 



P = 2 



dx 
v(x) 



(Al) 



where v(x) is the classical electron velocity at energy E x 



v(x) = \—{E x 



U(x)). 



(A2) 



The flux A is then given by the frequency of hits against 
the barriers, 2/P, times the transmission probability T 
across a barrier, and therefore: 



dx 
v 



(A3) 



This expression is very convenient because the transmis- 
sion coefficient eq. ^ for our parametrized potential, 
U(x), is known analytically |13|. For more general bar- 
rier profiles and the long lifetimes of interest, one can 
use the WKB approach and its corresponding connec- 
tion formulae across the barrier (see e.g. Appendix D of 
I): 



WKB 



k, dx = 



2rU tTTI \ 

-T (U(x) 



E x ) dx (A4) 



If the WKB wave function is used inside the well to de- 
termine t he p eriod P, the same decay half-life is obtained 
as in eq. |A3| above. 

Since the dot is located inside a crystal at temperature 
T", via phonon coupling the electrons in the dot should 
also be at the same temperature. The level occupations 
f(E) are determined by Fermi statistics: 



f(E) = [1 + e »bt> ] 

where these are now 2D energies, 
obtained from 



(A5) 

The Fermi level is 



(A6) 



i=(n a: ,ny) 



where the 2 accounts for spin degeneracy. For the ensem- 
ble of electrons in the dot, the flux A will now be the sum 
of fluxes for each occupied single particle level, weighted 
by the level occupancy: 



A = 



E V{E l )h 



(A7) 



i=(na: ,n y ) 



and the corresponding half-life is still r = 1/A. In par- 
ticular this argument applies in the T' = limit, as we 
implicitly assumed in Section II to explain the sequence 
of lifetimes. 



APPENDIX B: LIFETIME DEPENDENCE ON Q 

For a level of given n x , the lifetime depends on Q be- 
cause the barrier characteristics change, and so does the 
level energy, E nx . The latter varies mainly because Uq 
depends on Q, and this affects the transmission proba- 
bility T . To good approximation 



dE nx 
dQ 



dU 
dQ 



(Bl) 
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Neglecting the dependence of the level lifetime on the 
period P, we can write 



logic Tn. (Q) = log 10 r nsc (Q ) - O.U(Q - Q ) . (B5) 



d In 7 
~dQ~ 



rflnT dE 
dE dQ ' 



(B2) 



Taking the transmission probability T from the WKB 
expression leads to 



dlnr d 
dQ = dQ 



barri 



bar ri 



/ 2m* 



- E) dx 
( dU{x) 



h 2 U(x) -E\ dQ 



dE 
dQ 



dx 



(B3) 



Noting eq. Bl, the second contribution to the integral 
depends linearly on the placement of the potential floor. 
However, the variation of the barrier shape (dU(x)/dQ) 
cannot be neglected. Indeed, eq. @ gives approximately 



dUjx) 
dQ 



dU(x) 
dQ 




sinh 



e~ 2 ^ cosh 2 



when x < Xt, 




sinh 



3 e 2 ^ cosh 2 ( 
when x > Xb ■ 




Using |B4i the contribution from dU/dQ to the integral 
of eq. [B3| is obtained with an accuracy better than 2%. 




J. Cooper, C.G. Smith, D.A. Ritchie, E.H. Linfield, Y. 
Jin and H. Launois, " Direct Observation of Single Elec- 
tron Decay from an Artificial Nucleus", Proceedings of 
EP2DS-13, Physica E 6 (2000) 457-460. 
G. Gamow, Z. Phys. 51 (1928) 204. 

E.U. Condon and R. W. Gurney, Nature 122 (1928) 439. 
M.A. Preston and R.K. Bhaduri, "Structure of the Nu- 
cleus" Addison Wesley (1975), see Part III and Appendix 
D. 

W. van Dijk and Y. Nogami, Phys. Rev. Lett. 83 (1999) 
2867-71. 

J. Martorell and D.W.L. Sprung, to be published, and 

unpublished QUADRANT progress reports. 

J. Martorell and D.W.L. Sprung, Phys. Rev. B 49 (1994) 

13750-59. 

J. Martorell, Hua Wu and D.W.L. Sprung, Phys. Rev. B 
50 (1994) 17298-308. 

W. Monch, "Semiconductor Surfaces and Interfaces", 
Springer (1995). 

J. H. Davies, Semicond. Sci. Technol. 3 (1988) 995. 
J.H. Davies, A. Larkin and E.V. Sukhorukov, J. Appl. 
Phys. 77 (1995) 4504-12. 

J. Martorell and D.W.L. Sprung, Phys. Rev. B 54 (1996) 
11386-96. 

"Methods of Theoretical Physics", P.M. Morse and H. 
Feshbach, McGraw Hill, (1953); see pp. 1651-60. 
J Killingbeck, J. Phys. A 10 (1977) L99-L103; Phys. Lett. 
A 78 (1980) 235-6. 



FIG. 10. Dependence of lifetime on Q for the occupied 
levels with n x — 11: + signs; n x — 12: x signs; n x = 13: stars 
and n x — 14 open squares. The dash dotted lines correspond 



to eq. B5 



For the standard choice of parameters, and Q in the range 
300 to 310 the computed values of c?log 10 r/dQ turn out 
to be ~ —0.14 for the levels of interest. In Fig. [l0| we 
plot the evolution of the level lifetimes with Q, compared 
to the expression (Q = 303) 







